Low-threshold analysis of CDMS shallow-site data 
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Data taken during the final shallow-site run of the first tower of the Cryogenic Dark Matter Search 
(CDMS II) detectors have been reanalyzed with improved sensitivity to small energy depositions. 
Four ~224 g germanium and two ~105 g silicon detectors were operated at the Stanford Underground 
Facility (SUF) between December 2001 and June 2002, yielding 118 live days of raw exposure. Three 
of the germanium and both silicon detectors were analyzed with a new low-threshold technique, 
making it possible to lower the germanium and silicon analysis thresholds down to the actual trigger 
thresholds of ~lkeV and ~2keV, respectively. Limits on the spin- independent cross section for 
weakly interacting massive particles (WIMPs) to elastically scatter from nuclei based on these data 
exclude interesting parameter space for WIMPs with masses below 9GeV/c 2 . Under standard 
halo assumptions, these data partially exclude parameter space favored by interpretations of the 
DAMA/LIBRA and CoGeNT experiments' data as WIMP signals, and exclude new parameter 
space for WIMP masses between 3GeV/c 2 and 4GeV/c 2 . 
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I. INTRODUCTION 

Astrophysical evidence strongly suggests that matter 
constitutes approximately one quarter of the energy den- 
sity of the Universe. Baryons in stars and intergalactic 
gas account for only a small fraction of the matter den- 
sity, while the majority of the Universe's matter is of an 
unknown composition, collectively termed dark matter 
(DM) due to its apparently nonluminous nature [l| . Ob- 
servations of large-scale structure and supernovae, com- 
bined with measurements of the cosmic microwave back- 
ground, imply a total matter and energy density close or 
equal to the critical density Q. In terms of the critical 
density, the matter and energy budget breaks down as 
follows: 



Q haIyon = 0.0456 ± 0.0016, 



n UM = 0.227 ±0.014, 



and 



n 7OO+0.015 
U -' 28-0.016' 



(1) 



(2) 



(3) 



where fl\ represents the mysterious dark energy thought 
to be responsible for the current accelerating expansion 
of the Universe [H . 

The standard model of particle physics provides a sin- 
gle candidate for this nonbaryonic dark matter: the neu- 
trino. Large- volume neutrino observatories have success- 
fully measured and confirmed the existence of neutrino 
mass [4|. Nevertheless, neutrinos make only a small con- 
tribution to the dark matter density, 



93.14 eV/i 2 



< 0.006, 



(4) 



where m v < 0.28 eV and h is the dimensionless Hubble 
parameter [!, H|. Furthermore, formation of large-scale 
structure in the Universe constrains the neutrino com- 
ponent of the dark matter [f|. Simulations of structure 
formation require a significant nonrelativistic, or "cold," 
dark matter density |7j , which cannot arise from standard 
model neutrinos. 

A number of suitable dark matter candidates arise 
from theories that propose physics beyond the standard 
model. Weakly interacting massive particles (WIMPs) [8] 
are the most studied class of such dark matter parti- 
cles. In particular, many i?-parity-conserving weak-scale 
supersymmetric (SUSY) theories offer a natural dark 
matter candidate in the form of the lightest superpart- 
ner [l], 0, [l(|, often a neutralino. Massive, electrically 
neutral, and stable, the lightest neutralino, Xii of many 
SUSY theories is an excellent WIMP candidate. SUSY 
theories contain a vast space of unknown free parameters, 
which are constrained by requiring consistency with ex- 
isting empirical particle physics and astrophysics knowl- 
edge. Popular techniques for additionally restricting the 
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FIG. 1. (color online). Expected differential event rates for 
5 GeV/c 2 (top panel) and 100 GeV/c 2 (bottom panel) WIMPs 
scattering from Ge (blue/solid) and Si (red/dashed) targets. 
All event rate calculations are based on the standard halo 
model described in [3], for an arbitrarily chosen WIMP- 
nucleon cross section of lxlO -41 cm 2 . Each energy spectrum 
cuts off abruptly at a maximum recoil energy due to the as- 
sumed galactic escape velocity. We use the 544km/s galactic 
escape velocity from [l^], while all other halo parameters are 
taken from [lq |. and the local WIMP density is assumed to 
be 0.3GeV/cm 3 . 



extent of the free parameter space result in a lower bound 
on the Xi mass of ^40GeV/c 2 [ll|. Efforts to explore 
a wide range of SUSY free parameter space indicate 
that Xi masses as low as a few GeV/c 2 can be accom- 
modated [l2| . Under some scenarios, a relatively light 
WIMP could resolve the apparent conflict between the 
DAMA/LIBRA annual modulation signal and the null 
results of other experiments 131 ]. 

If WIMPs are the dark matter, they form a spherical 
cloud (or halo) in which the luminous portions of the 
Milky Way are embedded, and will scatter very rarely off 
the nuclei in terrestrial matter. Direct detection experi- 
ments seek to observe and measure the kinetic energy, in 
the keV range, of the recoiling nuclei. The expected spec- 
trum of WIMP-induced nuclear recoils decreases rapidly 
with increasing recoil energy, with the mean recoil en- 
ergy directly proportional to the reduced mass of the 
WIMP-nucleus system. Events with the lowest recoil en- 
ergies are most numerous for all WIMP masses, and thus 
direct detection experiments generally strive for low re- 
coil energy thresholds. Sensitivity to low recoil energies 
is particularly crucial for experiments seeking to detect 
light WIMPs. Figure [1] illustrates expected spectra of 
nuclear-recoil energies for two types of target nuclei for 
two WIMP masses. 

The advantages of a low threshold must be weighed 
against the degradation of background rejection capabil- 
ity at low recoil energies. Furthermore, the intrinsic rates 
of some categories of background may also increase at low 
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energies. Finally, special care is required to account for 
the effects of non-zero energy resolution at energies near 
the electronic-noise level. 

In this work we describe a new analysis of data from 
the Cryogenic Dark Matter Search (CDMS) experiment, 
with special attention to events of low recoil energy. In 
order to access this low-energy parameter space, we forgo 
the pulse-shape discrimination techniques used to reject 
near-surface background events in previous CDMS anal- 
yses. In contrast to previous CDMS results, the signal 
region of this analysis will be populated by a number 
of background events. We adopt an inclusive philoso- 
phy that maximizes the detection efficiency at low en- 
ergy while limiting the rate of nonphysical sources of 
background, such as electronic noise. This low-threshold 
analysis sacrifices some of the strengths of the CDMS ex- 
periment's traditional background discrimination meth- 
ods for a chance to probe previously untested low-mass 
WIMP parameter space. 



II. EXPERIMENTAL SETUP AND TECHNIQUE 
A. Apparatus 

The data reported on here were recorded during the 
final exposure of the first six CDMS II detectors at the 
Stanford Underground Facility (SUF) [U The SUF 
setup provided a high level of shielding against external 
sources of radiation. The SUF was a shallow site with 
~17 meters water equivalent overburden, effectively stop- 
ping hadronic cosmic rays and reducing the muon flux by 
a factor of 5. The remaining incident muons were tagged 
with a high-efficiency, hermetic plastic scintillator muon 
veto, allowing offline rejection of muon-coincident detec- 
tor interactions. The muon veto enclosed several layers 
of tightly packed passive shielding. A 15 cm-thick outer 
lead shield and 25 cm-thick outer polyethylene shield sur- 
rounded the detector cold volume to attenuate external 
photons and degrade external neutrons, respectively. In- 
side the radiopure copper walls that delineated the in- 
nermost 20 mK cold volume, 1 cm of ancient lead and an 
additional 11kg of polyethylene surrounded the detector 
assembly, providing further shielding. 

In the center of the apparatus six Z-sensitive 
Ionization- and Phonon-mediated (ZIP) detectors [l9j 
were arranged in a vertical stack ("tower"), with adja- 
cent detectors separated by 2.2 mm with no intervening 
material. Table|T]indicates the detector names, materials, 
masses and relative positions within the tower. 

Each detector has two ionization electrodes deposited 
on its bottom surface; a circular inner electrode ("q- 
inner") covering 85% of the physical area, and an an- 
nular outer electrode ("q-outer") which permits identifi- 
cation and rejection of events with energy depositions 
within the outer detector volume. Photolithographed 
onto the top side of each Ge (Si) detector were 4144 
(3552) Al and W superconducting Quasiparticle-trap- 



TABLE I. The first six CDMS II detectors are listed in order 
of their relative positions (from top to bottom) within the 
detector tower, indicating each detector's name, material, and 
mass. 



Name Material Mass (grams) 

Zl Ge 230.5 

Z2 Ge 227.6 

Z3 Ge 219.3 

Z4 Si 104.6 

Z5 Ge 219.3 

Z6 Si 104.6 



assisted Electrothermal-feedback Transition-edge sensors 
(QETs). The 1036 (888) QETs in a given Ge (Si) detec- 
tor quadrant are electrically connected, resulting in four 
individually read out phonon sensors whose shared bor- 
ders orthogonally bisect the surface. By measuring tim- 
ing and pulse height differences between the sensors, we 
can reconstruct an event's position in the plane parallel to 
the detector's top and bottom surfaces ( "xy-position" ) . 

Following amplification by a SQUID array [2 6] and 
room temperature electronics, two copies of each detec- 
tor's four phonon signals were generated at the hard- 
ware level. A band-pass filtered analog sum of one set 
of phonon signals (the "triggering phonon energy") was 
compared to a low-level discriminator threshold. The 
resulting logical pulses were OR'ed across all six detec- 
tors to form the experimental trigger. The second set of 
phonon signals were individually digitized, subjected to 
a software optimal filter, and summed to constitute the 
"reconstructed phonon energy." The hardware band-pass 
filter's poles were chosen to resemble the software filter- 
ing as closely as possible. The energy deposition used to 
trigger the data acquisition system is therefore very sim- 
ilar to the energies evaluated in software during offline 
analysis, but not exactly the same. The discriminator 
thresholds were carefully tuned to the lowest levels possi- 
ble to ensure that the overall trigger rate of <1 Hz during 
WIMP search was dominated by true particle interac- 
tions, while simultaneously allowing occasional triggers 
due to electronic-noise fluctuations. 



B. Measurement technique 

Energy deposited by recoils causes two types of sig- 
nals in our detectors. Most of the energy is deposited 
as a spectrum of high frequency athermal phonons. In 
addition, electron-hole pairs arc created. The average de- 
posited energy per pair created is e=3.0 (3.8) eV for Ge 
(Si) [2l| (see also Appendix C in [22| for a detailed dis- 
cussion of e). To measure the ionization, the QET side of 
our detectors was held at ground, and a bias voltage be- 
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tween -3 V and -6 V was applied to the electrodes on the 
other detector side, causing drift of liberated electrons 
(holes) to the QET (electrode) side. Charge-sensitive 
amplifiers with field-effect transistor front-ends amplified 
the signals on the electrodes and enabled inference of 
the "ionization energy" Q from the total liberated ion- 
ization. The total phonon signal is comprised of three 
parts: l)phonons produced promptly by the recoil ("pri- 
mary phonons"); 2)phonons produced by the driftin g io n- 
ization ( "Neganov-Luke phonons" or "drift heat" [23jj ) ; 
and 3)phonons produced when electrons and holes re- 
combine, usually at the detector surfaces ("recombina- 
tion phonons"). The QETs and their amplification sys- 
tem enabled deduction of the energy Ptotai from the total 
phonons. The "recoil energy," £? r ecoil> is equal to the sum 
of the primary and recombination phonon energies, and 
can be deduced by subtracting off the drift heat: 

eV 

Precoil = Ptotai Q, (5) 

e 

where V is the absolute value of the bias voltage, and e 
is the charge of the proton. 

The CDMS detectors use these two sensor technolo- 
gies to discriminate nuclear recoils, produced by WIMP 
candidates, from much more numerous electron recoils, 
produced mostly by background photons. The detectors 
provide a simultaneous measurement of ionization and 
phonons for each particle-interaction event within the 
target. The ratio of ionization to recoil energy ("ion- 
ization yield" Y = Q/E vccol \) is higher for electron re- 
coils than for nuclear recoils, and provides near-perfect 
event-by-event discrimination for recoil energies in ex- 
cess of lOkeV. However, the ionization yield is broadened 
by electronic noise for recoil energies lower than lOkeV, 
and discrimination power is lost for recoil energies below 
about 2keV. 

The energy scale for Q and P roC oii was calibrated by 
electron recoils caused by gamma-ray sources. The re- 
sulting ionization yield for electron recoils is (on aver- 
age) equal to unity by construction. However, nuclear re- 
coils cause less ionization, and generally follow the theory 
by Lindhard et al. [241 ] ■ Nuclear recoils from a neutron 
source allowed determination of the energy dependence 
of their ionization yield, which smoothly increases from 
-0.2 at E iecoi i=2keV to -0.4 at 100 keV. 

The ionization yield for nuclear recoils with recoil en- 
ergies less than 2 keV is difficult to measure directly with 
ZIP detectors. However, measurements using Ge detec- 
tors have been performed for ener gies as low as ~0.25 keV 
(see Fig. 3 in Appendix III in [25] for example), and 
agree well with Lindhard et al.'s model. We therefore 
extrapolate via a power law from 2 keV to lower ener- 
gies such that the ionization yield is zero at keV. The 
extrapolation is needed only to estimate the energy de- 
pendence of a few detection efficiencies for the lowest- 
energy nuclear recoils. At these low energies, the recoil 
energy is less than (but nearly equal to) Ptotai, and the 
extrapolation is used to estimate recoil energy from Ptotai 



via a small correction. Any uncertainty related to this 
small correction is therefore small as well. Alternative 
extrapolation methods have been used to test the effect 
that uncertainty in the ionization yield for low-energy 
nuclear recoils might have on our ability to detect light 
WIMPs. We have determined that any systematic uncer- 
tainty introduced by our particular choice is small com- 
pared to the other sources of uncertainty discussed in 
Section |IVDj 

As discussed earlier, in this paper we extend the lowest 
recoil energies to the IkeV (2keV) range for Ge (Si) tar- 
gets. Nuclear recoils with such low energies will appear to 
deposit very small ionization energy, in the tenths of keV. 
Although earlier CDMS analyses implemented minimum 
requirements on the ionization energy, in this analysis 
we have avoided any minimum requirement on Q. We 
emphasize that our hardware trigger used only the total 
phonon energy, which slightly exceeds the recoil energy 
for nuclear recoils. 



C. Recoil-energy estimators 

One way to estimate the recoil energy is with both 
measured Ptotai and measured Q used in Eq. [5] on an 
event-by-event basis. This method incorporates errors 
on both these measured quantities into the estimate of 
Procoib and introduces significant correlations between 
-EVecoii and the ionization yield. This is the traditional 
CDMS method for estimating recoil energy because it is 
accurate even if it is not known whether an event is an 
electron or nuclear recoil. The analysis cuts and their 
efficiencies (described in Section IlII B[) use this event- by- 
event recoil-energy estimate when needed. We refer to 
this estimate as the "Q-corrected" recoil energy to dis- 
tinguish it from an alternative method discussed below. 
Throughout this paper, the recoil energy is always Q- 
corrected unless otherwise stated. 

The recoil energy can also be estimated from Ptotai 
alone by scaling it to reflect the average ionization yield 
response measured from calibration samples. For exam- 
ple, since the ionization yield for electron recoils is on 
average equal to one, the ionization energy for electron 
recoils is on average equal to the recoil energy itself. Re- 
placing Q by Procoii in Eq. [5] yields P re coii equal to one- 
half (one-third) of Ptotai for Ge ZIP detectors operated 
with a -3 V (-6 V) bias voltage. We refer to this estimate 
as the recoil energy corrected by electron-recoil ioniza- 
tion yield, or ">ER-corrected" recoil energy. Due to its 
superior spectral resolution, it is particularly useful when 
studying ZIP-detector response to x-ray and gamma-ray 
sources (see Fig. [6] for example). 

A similar recoil-energy estimate can be made for nu- 
clear recoils. The reduced drift heat associated with 
nuclear recoils is energy dependent, and can be sub- 
tracted according to the mean ionization yield measured 
from 252 Cf calibrations (dashed line in Fig. |4] for ex- 
ample). The resulting recoil-energy estimator has the 
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advantage of including electronic noise from only the 
phonon channel. We refer to this estimate as the re- 
coil energy corrected by nuclear-recoil ionization yield, 
or "iNR-corrected" recoil energy. Discussed in more de- 
tail in Section IIII A[ the hardware and software energy 
thresholds employed in this analysis depend solely on the 
total phonon signal. The associated threshold efficien- 
cies (see Figs. [2] and El and Table [TT] for example) for 
detecting nuclear recoils are therefore functions of the 
iNR-corrected recoil energy. 



D. Data samples 

The first tower of CDMS II ZIPs was commissioned at 
the SUF in the second half of 2001, and WIMP-search 
data were recorded between December 2001 and June 
2002. Collectively termed "Run 21," the WIMP search 
was split into two distinct operational periods. The Ge 
(Si) detectors were initially operated with a -3V (-4V) 
bias voltage for 94 days between December 2001 and 
April 2002, yielding just over 66 live days (henceforth 
referred to as the "3V data"). These data were origi- 
nally analyzed with a 5 keV analysis threshold, and the 
resulting exclusion limit on the spin-independent WIMP- 
nucleon cross section is still one of the strongest con- 
straints on low- mass WIMPs [l?} ■ An additional 52 live 
days of previously unpublished data were recorded over a 
period of 74 days between April and June 2002. In an at- 
tempt to reduce the surface-event background, we exper- 
imented with a -6V bias voltage during this latter part 
of the run (henceforth referred to as the "6V data" ) . Al- 
though the larger bias voltage improved charge collection 
for events near the detector surfaces, it also degraded the 
phonon pulse rise-time information used to reject surface 
events. Overall, the rate of surface events leaking into the 
signal region following a surface-event rejection analysis 
was slightly greater for the 6V data. CDMS ZIP detec- 
tors have since been operated with the lower bias voltage 
settings. In this paper we analyze data from both charge 
bias runs. 

A mixture of WIMP-search and calibration data was 
recorded during Run 21. Detector response to elec- 
tron recoils was tested by introducing a 60 Co gamma- 
ray source just inside the outer lead shielding. Three 
gamma-ray calibrations were performed for each volt- 
age bias. Relatively short exposures conducted before 
recording any WIMP-search data were followed by more 
extensive exposures midway through and after WIMP- 
search runs. These data were primarily used for an event- 
by-event correction of the phonon signals for observed 
position- and energy- dependent response functions ( "po- 
sition correction"). The xy-position information avail- 
able from the modularity of the phonon sensors was used 
to correct variations in the reconstructed optimal-filter 
pulse height due to variation in pulse shape with event 
position, yielding improved phonon energy resolution. 

With part of the outer polyethylene shielding removed, 



the detectors were also exposed to a 252 Cf neutron source 
to calibrate their response to nuclear recoils. Neutron 
calibrations were performed at the beginning and end of 
the WIMP searches for each bias voltage. The ioniza- 
tion yield distributions obtained from these data helped 
define the nuclear-recoil signal region that WIMPs are 
expected to populate. Nuclear-recoil detection efficien- 
cies were estimated from these data as well. 

Although the hardware thresholds were tuned to oc- 
casionally trigger on electronic noise, the rate was far 
too low for proper characterization of the distribution 
of subthreshold noise pulses, and for estimates of near- 
threshold efficiencies. To understand the near-threshold 
phenomena, samples of detector performance without the 
trigger bias ( "nontriggered data") are critically impor- 
tant. During Run 21, nontriggered data were obtained 
by invoking a randomly generated software trigger, re- 
sulting in a periodic sampling of the electronic-noise en- 
vironment. These events constitute ~5% of all triggers, 
and can be considered a third type of calibration data. 
Due to the low rate of true particle interactions in the 
detectors, randomly triggered events usually consisted 
of only noise fluctuations, with a very low probability 
of measuring a non-zero energy in any of the detectors' 
sensors. While useful for determining daily noise levels, 
the randomly triggered data are inadequate for efficiency 
estimates since so few have energies near the detector 
thresholds. Furthermore, due to a timing artifact asso- 
ciated with the implementation of the random trigger, 
these data were plagued by a low rate of unrepresenta- 
tive noise traces. When reconstructed, the problematic 
traces contribute non-Gaussian tails to each detector's 
underlying electronic-noise distribution. Consequently, 
if the randomly triggered noise distributions are scaled 
to an exposure equivalent to the WIMP-search data and 
then subjected to the WIMP-search cuts, the resulting 
rate of noise pulses in the signal region is grossly overes- 
timated. We therefore made no use of nontriggered data 
for characterizing resolutions and efficiencies. 

Fortunately, an alternate sample of data that avoids 
the trigger bias was acquired during normal, triggered 
readout. Although most triggers were associated with 
an energy deposition in a single detector, traces for all 
six detectors were recorded and analyzed. Data from the 
detectors that did not trigger ( "other-detector triggers" ) 
provide a fair and representative sample of each detec- 
tor's subthreshold noise distribution. Provided the noise 
distribution is constructed from events where the trigger- 
ing detector was not an adjacent detector, non-Gaussian 
tails are avoided. Events for which the triggering detector 
was either directly above or below often contained small 
energy depositions due to true multiple-detector interac- 
tions. This class of events includes a sampling of recoil 
energies up to and exceeding the detector trigger thresh- 
olds, providing events with which to probe near-threshold 
behavior. Some of these had energies exceeding their de- 
tector's trigger threshold and are not truly nontriggered 
events, but rather events with a delayed trigger that oc- 
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FIG. 2. (color online). Hardware trigger efficiency estimate 
(error bars) for a representative Ge detector (Z2 3V data) 
as a function of YNR-corrected recoil energy. The efficiency 
is calculated in bins of 0.25 keV by dividing the distribution 
of energies correlated to the presence of logical pulses in the 
post-trigger history (green/light solid) by the distribution of 
all reconstructed energies (black/dark solid). The efficiency 
scale is given by the y-axis on the right (blue/dash-dotted 
grid lines), while the scale for the histograms is given by the 
y-axis on the left (black/dotted grid lines). A split- width error 
function (red/dashed) fits the efficiency estimate over the full 
range of energies, yielding a hardware threshold of 0.74 keV 
at 50% efficiency. For this detector the hardware trigger is 
100% efficient for energies above ~3keV. 



curred during a readout instigated by another detector. 
A logical-pulse-based post-trigger history of when these 
delayed triggers occurred was recorded for each detec- 
tor and event, and is the basis of the hardware trigger 
efficiency estimates detailed in the next section. 



III. DATA ANALYSIS 

A. Noise and thresholds 

To be considered viable WIMP candidates, events were 
required to exceed hardware and software thresholds in 
phonon energy. We used data from other-detector trig- 
gers to measure the corresponding detection efficiencies 
as a function of reconstructed phonon energy. Since the 
triggering phonon energy was slightly different from its 
offline reconstructed counterpart, the hardware thresh- 
old efficiencies as a function of the latter are not simple 
step functions. To characterize the hardware threshold 
efficiencies, we first evaluate the reconstructed phonon 
energy in the 50 /is following the time of each other- 
detector trigger. Generally these reconstructed phonon 
energies are Gaussian distributions centered at zero en- 
ergy, consistent with electronic noise ( "noise cores" ) . Oc- 



TABLE II. The phonon energy thresholds (at 50% efficiency) 
with la errors (68.3% C.L.) for the viable low-threshold detec- 
tors in terms of YNR-corrected recoil energy. For each detector 
the hardware trigger threshold and the average 6er noise re- 
jection software threshold is listed for its 3V and 6V data. 



Detector 


Hardware (keV) 
3V 6V 


Software (keV) 
3V 6V 


Z2 
Z3 
Z4 
Z5 
Z6 


n 74+0-07 

, -iq+0.07 

1 yy+0.12 
L -< 1 -0.10 

1 00+ 05 

1 cq+0.08 
x - JO — 0.06 


67+0- 09 

U - D ' -0.02 

-i 19 +0.07 
L - lz -0.05 

-1 71 +0.13 
L -' x -0.10 

qi+0-06 

1 cc+0.09 
l-JJ-0.06 


0.63 ±0.01 

0.82 ±0.02 
1 fi2 +o.io 

0.73 ± 0.02 

1 -JQ+0.06 
1 - oa -0.05 


0.58 ±0.02 
0.72 ± 0.03 

i t-o+0.09 
1 - oz —0.0B 

fi Q+0.04 
u - oo -0.03 

1.34 ±0.05 



casionally nonzero energies were reconstructed, augment- 
ing each noise core with an approximately uniform dis- 
tribution extending to higher energies. Sometimes the 
presence of a logical pulse from the low-level discrimina- 
tor threshold was recorded in the post-trigger history of 
an other-detector trigger. We used events with logical- 
pulse presence to obtain the distribution in reconstructed 
phonon energy of successful triggers, and then divided by 
the distribution of all (unsuccessful and successful) trig- 
gers to characterize the hardware threshold efficiency, as 
illustrated for one Ge detector in Fig. [2] 

A reconstructed phonon energy threshold ("software 
threshold" ) was applied based upon the essentially Gaus- 
sian behavior of the reconstructed phonon energy. A 
Gaussian fit to the noise core of each day's randomly trig- 
gered events was performed, and we required that events 
exceed a threshold of 6a above the mean (for that day) in 
reconstructed phonon energy to be considered a WIMP 
candidate. The 6a thresholds for the Z2 and Z3 3V data 
were the most stable, varying by less than 5% from day 
to day, while those for the Zl data were the least sta- 
ble, varying by as much as 50%. The software threshold 
efficiencies were calculated using the same method em- 
ployed for the hardware trigger efficiency estimates de- 
scribed above, and differ from step functions because of 
time variations in the Gaussian fit parameters. Table |TT] 
summarizes the average hardware and software phonon 
energy thresholds. 

The topmost Ge detector (Zl) is not included in Ta- 
ble pi] because it was rejected as a low-threshold detector. 
Zl suffered from particularly strong dependence of pulse 
height on xy-position for which the gamma calibration 
position correction was unable to fully compensate. The 
corrected energy resolution and associated 6a thresholds 
are three to four times larger than those of the other Ge 
detectors, yielding not only a larger analysis threshold, 
but also causing the ionization yield-based discrimina- 
tion to break down at a higher energy. Since a detector's 
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FIG. 3. (color online). Fits to binned estimates of the hard- 
ware (orange/light solid) and software (blue/dashed) phonon 
energy threshold efficiencies are multiplied, yielding the com- 
bined threshold efficiencies (black/dark solid) for typical Ge 
(top panel, Z5 6V data) and Si (bottom panel, Z4 3V data) 
detectors as a function of iNR-corrected recoil energy. 

low-mass WIMP sensitivity is critically dependent on its 
detection threshold, data from Zl would contribute very 
little to the reach of this analysis while allowing a dispro- 
portionate number of background events to leak into the 
signal region. Zl was only used to veto multiply scatter- 
ing events. 

Each detector's final threshold efficiency is a product 
of its hardware and software threshold efficiencies. As 
a function of reconstructed phonon energy, the software 
threshold efficiencies nearly resemble step functions, ris- 
ing to 100% more quickly and at lower energies than the 
hardware efficiencies. Consequently, the combined effi- 
ciencies differ from the hardware trigger efficiencies only 
for very low energies. Figure [3] illustrates the combined 
efficiency for typical Ge and Si detectors. For all but one 
detector the combined thresholds at 50% efficiency are 
equal to the hardware thresholds listed in Table [TTJ Dur- 
ing one week of the 3V data run, Z4's phonon noise was 
abnormally high, causing an extended period of higher 
software thresholds. The effect on the software and com- 
bined efficiencies is shown in Fig. [3j The resulting com- 
bined threshold of 1.9±0.1keV is the largest among the 
accepted low-threshold detectors. 



B. Analysis cuts and efficiencies 

WIMP candidates were required to pass a variety of 
stringent data-quality cuts. Events for which any part 
of the data record was incomplete or inconsistent were 
excluded. Excessively noisy traces as well as traces with 
multiple pulses ( "pileup" ) were rejected based on the per- 
formance of the optimal-filter fits to the ionization sig- 



nals. Fits with unusually large x 2 values were removed 
from consideration. Noisy traces and pileup were further 
suppressed by requiring the variances of the phonon and 
ionization traces' prepulse baselines to be within 5a of 
the average behavior for randomly triggered noise traces. 
Together, the data-quality cuts exhibit no energy depen- 
dence, and have a combined WIMP detection efficiency 
of ^99%. This applies for each detector type and for 
both the 3V and 6V data. 

If this analysis were restricted to recoil energies greater 
than 2keV, the above data-quality cuts would be suffi- 
cient. Throughout the 6V data run, however, the sec- 
ond Ge detector (Z2) was subject to intermittent peri- 
ods of high trigger rates ("event bursts") due to sub- 
2 keV pulses. Many followed cryogenic and detector ser- 
vicing periods during which detector temperatures were 
elevated. Consequently, we conclude that these event 
bursts were not due to physical recoils in the detector, 
and were considered periods of poor data quality. Close 
examination of the phonon traces for these events re- 
vealed elevated prepulse baselines with nominal variabil- 
ity, allowing them to pass the standard data-quality cut 
described above. Discrimination parameters based on 
the traces' average prepulse baselines were developed to 
cut the event bursts. Application of the cut reduces Z2's 
6V WIMP search from 52 days to 20 days, while reduc- 
ing the corresponding number of WIMP candidates by a 
factor of ~30. Despite being isolated to a single detec- 
tor and less than half of the WIMP-search data, without 
this cut the burst events would be the dominant source 
of background invading the signal region. 

A fiducial-volume cut based on the ionization signals 
rejects recoils that occurred near detector edges. The q- 
inner signal estimates an event's ionization energy, while 
the q-outer signal is required to be consistent with noise. 
Crosstalk between the inner and outer electrodes adds a 
small contribution to the q-outer signal proportional to 
the q-inner pulse amplitude. To model this dependence, 
the q-outer noise levels were parametrized as a function 
of q-inner energy, yielding 2a q-outer "noise bands" de- 
signed to accept 95.5% of events within the fiducial vol- 
ume. Phonon sensors of the adjacent detector also in- 
duced signals in the electrodes, initially causing a loss 
of otherwise viable events. To retain these events, the 
fiducial-volume cut was modified to reject only events 
that exceed their noise band's upper limit, increasing the 
cut's expected acceptance to 97.7%. Combining the ac- 
ceptance with the physical coverage of the q-inner elec- 
trode yields an expected fiducial volume of ^83% of the 
gross detector mass. We measured the efficiency of the 
fiducial- volume cut for recoil energies between 4 keV and 
100 keV with nuclear-recoil calibration data from 252 Cf 
exposures. The average efficiency agrees well with physi- 
cal expectation, varying between ~81% (Z4 3V data) and 
^83% (Z2 3V data). Over the measurable energy range, 
the efficiency of the cut exhibits a weak energy depen- 
dence, generally decreasing with increasing recoil energy. 
The Z4 6V data are the most extreme, where the effi- 
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ciency smoothly decreases from ^84% at 4keV to ~79% 
at 100 keV. For recoil energies less than 4keV, a nuclear- 
recoil event's q-inner and q-outer signals are difficult to 
distinguish from electronic noise, and the fiducial- volume 
cut is unable to differentiate events within the fiducial 
volume from those that occurred near a detector's outer 
edge. Due to increased acceptance, the efficiency of the 
fiducial-volume cut should therefore increase rapidly as 
the recoil energy decreases toward keV. Since we were 
unable to measure this low-energy behavior, we made the 
conservative choice to linearly extrapolate the efficiency 
for recoil energies less than 4 keV to match the efficiency 
measured at 4keV. 

The energy deposited by a WIMP in one of our detec- 
tors would be so localized and infrequent that only events 
caused by backgrounds will cause significant energy de- 
positions in two or more detectors simultaneously. We 
therefore impose a "single-scatter" criterion, requiring 
signal events to have had a significant energy deposition 
in no more than one detector. Since the phonon signal 
provides the most sensitive indicator for a particle inter- 
action in our detectors, we based the single-scatter cut 
on the 6a software thresholds described above and in Ta- 
ble HU An event is considered a single scatter if only one 
detector has a reconstructed phonon energy exceeding its 
software threshold. With an experimental trigger rate of 
<lHz, the probability of more than one pulse occurring 
within our <2 ms digitization time is negligible. The effi- 
ciency of the single-scatter cut is therefore nearly 100%, 
with the near-threshold behavior taken into account by 
the software threshold efficiency described above and in 
Fig.H 

Candidate events must occur when there is no activ- 
ity in the muon veto. A time history of muon veto ac- 
tivity ("veto hits") was recorded for each ZIP-triggered 
event at an average rate of ~5 kHz. This rate far exceeds 
the true rate of incident muons. Due to large and awk- 
ward counter geometries, many of the plastic scintillator 
counters had regions with poor light collection. To en- 
sure efficient tagging of muons passing through these re- 
gions, the photomultiplier tubes monitoring each counter 
were operated at very high gains, causing regions with 
superior light collection (particularly near the photomul- 
tiplier tubes) to be partially sensitive to environmental 
gamma radiation. We achieved a muon tagging efficiency 
of >99.9% by rejecting a relatively high rate of ZIP de- 
tector interactions that were accidentally coincident with 
gamma rays registering as veto hits. An event is vetoed 
if there were any veto hits in the 50 /is to 80 /zs preced- 
ing the triggering phonon pulse, allowing for a significant 
difference in veto and phonon signal arrival times. The 
veto signals were effectively instantaneous relative to the 
more slowly rising phonon pulses. The longer delay was 
assigned to events with recoil energies <3.5keV, allow- 
ing smaller phonon pulses sufficient time to rise past the 
hardware trigger thresholds. WIMP detection efficiency 
is lost due to veto hits from gamma-ray activity. The de- 
tection efficiency that remains for each detector following 




Recoil energy (keV) 



FIG. 4. (color online). Ionization yield plotted as a func- 
tion of recoil energy for representative Ge (Z3 6V data) 252 Cf 
neutron calibration data, with the average ionization yield 
(blue/dashed), 2a nuclear- recoil band (blue/dark solid) and 
software phonon threshold (yellow/light solid) overlaid. Un- 
vetoed single scatters passing the fiducial-volume and data- 
quality cuts are displayed for events consistent with the 
nuclear-recoil criterion (orange/light dots), and some that are 
not (black/dark dots). A substantial gamma-ray flux from the 
252 Cf source populates the "electron-recoil band" near an ion- 
ization yield of 1. Two lines at ~10.4keV and ~66.7keV re- 
sulting from decays of Ge isotopes can be distinguished among 
the electron-recoil events. As described in Section [III CI the 
low number of events above the nuclear-recoil band between 
2 keV and 6 keV recoil energy qualitatively demonstrates that 
channeling does not significantly diminish this experiment's 
efficiency for detecting low-mass WIMPs. 



application of the veto cut ranges from ^67% to ^78%. 

The final selection criterion requires a candidate event 
to have an ionization yield consistent with being a nu- 
clear recoil, thereby discriminating against the otherwise 
overwhelming rate of electron recoils. Several short ex- 
posures to the fission neutrons from a 252 Cf source pro- 
vided sufficient nuclear-recoil events to parametrize each 
detector's ionization yield response as a function of re- 
coil energy. When binned according to recoil energy, a 
detector's ionization yield distribution for nuclear recoils 
is well described by a Gaussian for energies as low as 
2keV. The mean ionization yield and associated la width 
as a function of recoil energy were estimated in a binwise 
fashion with maximum likelihood fits. A simplified form 
of Lindhard et al.'s theory for nuclear-recoil ionization 
yield in semiconductor crystals (Y = a * -E^ecoU , where 
a and b are the fit parameters) was fit to the resulting 
collection of means, while an inverse-squared form was 
fit to the squared widths (ay = d 2 + c 2 /E 2 ccoil , where c 
and d are the fit parameters). The fitted ionization yield 
trends were used to construct nuclear-recoil acceptance 
bands, where an event is considered a nuclear recoil if 
its ionization yield lies within 2a of the mean. An ex- 
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FIG. 5. (color online). The detection efficiencies associated 
with the hardware and software phonon thresholds (top panel, 
Si Z6 3V data) and the combined analysis cuts (bottom panel, 
Ge Z2 6V data) for representative detectors. The best-fit 
mean efficiency (orange/light solid) is compared in each case 
to its 90% statistical lower limit (blue/dark solid). The com- 
bined threshold efficiencies are plotted versus iNR-corrected 
recoil energy, while the combined efficiency of the analysis 
cuts is plotted as a function of Q-corrected recoil energy. The 
step in the latter is due to the larger muon-veto window cho- 
sen for the low-energy events. Due to the uncertainty intro- 
duced by the efficiencies' statistical errors, we conservatively 
used the 90% lower limit (1.28ct below the mean) when cal- 
culating upper limits on the WIMP-nucleon cross section. 



ample nuclear-recoil band and mean ionization yield fit 
are compared to the 252 Cf events from which they were 
formed in Fig. @] The selection efficiency for each detec- 
tor's band is well fit by an energy independent efficiency 
ranging from -93% (Z6 3V data) to -96% (Z3 and Z4 
3V data). 

The combined efficiency of the analysis cuts as a func- 
tion of recoil energy for a representative detector is shown 
in Fig. [5] The step between 3 keV and 4 keV is due to 
the larger muon-veto window chosen for the low-energy 
events. The 90% confidence level (statistical) lower limit 
efficiency (1.28c below the mean) is also shown, providing 
an indication of the statistical accuracy of the efficiency 
estimate for the analysis cuts. 



C. Energy scale and resolution 

Our ZIP detectors are characterized by two energy 
scales: recoil energy for electron recoils, and recoil en- 
ergy for nuclear recoils. The electron-recoil energy scale 
was initially calibrated with 662 keV gamma rays from a 
137 Cs source as part of a series of short diagnostic and de- 
bugging runs conducted at the very beginning of Run 21. 
The gamma rays emitted by 137 Cs are lower in energy 
relative to the >1 MeV gamma rays emitted by 60 Co, 




Recoil energy (keV) [Y„ -corrected] 



FIG. 6. Event rate of electron recoils for representative Ge (Z3 
6V) WIMP-search data, plotted versus YER-corrected recoil 
energy. Included events are unvetoed single scatters passing 
the data-quality and fiducial-volume cuts. Three peaks due 
to internal decays of unstable Ge isotopes are clearly visible 
at 1.3 keV, 10.4 keV and 66.7 keV. 



allowing easier calibration of the Si detectors for which 
such energetic photons are only partially contained. 

Each Ge detector's electron-recoil energy scale was 
confirmed and monitored using three distinct energy 
peaks that conveniently span our analysis energy range. 
Both 68 Ge and 71 Ge are unstable isotopes; 68 Ge is pro- 
duced by cosmic rays, while 71 Ge by thermal neutron 
capture. The former is long lived with a half-life of —271 
days and was primarily produced during detector con- 
struction and testing at sea level, while the latter is short 
lived with a half-life of 11.43 days and was primarily pro- 
duced during exposure to the 252 Cf source [2y|. Both iso- 
topes typically decay through electron capture followed 
by emission of x-rays or Auger-electrons, of 10.4 keV for 
K-capture (-90%) or 1.3keV for L-capture (-10%). We 
measured the ratio of L- to K-captures to be 0.122±0.009, 
in good agreement with the ratio of 0.12 listed in [26| . 
The cosmogenically- induced meta-stable 73m Ge state de- 
cays through successive emission of 53.4 keV and 13.3 keV 
photons, resulting in a 66.7 keV peak at the upper end 
of our energy range. The spectrum for an example Ge 
detector shown in Fig. [6] clearly shows the three peaks 
from electron recoils induced by the above processes. To 
improve spectral resolution, Fig. [5] is plotted as a func- 
tion of YER-corrected recoil energy. A similar spectrum 
with corresponding features can be constructed from the 
ionization signal as well. 

The Ge lines allowed measurement of the Ge phonon 
and ionization energy resolutions across our full analy- 
sis energy range. The fractional energy resolutions for 
electron recoils vary from 8.8% (Z5 6V phonon chan- 
nel) to 26.0% (Z3 3V ionization channel) at 1 keV, and 
from 2.0% (Z3 and Z5 3V ionization channels) to 5.4% 
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TABLE III. The la recoil-energy resolution at zero energy 
for the viable low-threshold detectors is listed for both the 
3V and 6V data. The YNR-corrected resolution is based on 
the total phonon signal scaled to an equivalent recoil energy 
for nuclear recoils, while the Q-corrected resolution includes 
electronic-noise contributions from both the ionization and 
phonon signals. All resolutions are rounded to the nearest eV 
and have an accuracy of ~1%. 





Y NR - 


corrected (eV) 


Q-corrected (eV) 


Detector 


3V 


6V 


3V 6V 


Z2 


95 


88 


223 387 


Z3 


126 


114 


282 477 


Z4 


208 


196 


398 538 


Z5 


113 


102 


258 450 


Z6 


185 


179 


434 678 



(Z2 3V phonon channel) at 100 keV. As there are no 
peaked features due to internal Si isotope decays, the 
same could not be done for the Si detectors except at 
zero energy, where resolution was determined solely by 
electronic noise. 

A low-threshold analysis must take into account the 
effect of non-zero energy resolution near the detection 
threshold. The dominant contribution to the resolu- 
tion near threshold for all detectors was from electronic 
noise. We therefore dropped the energy-dependent terms 
of each detector's energy resolution when calculating 
WIMP exclusion limits, and simply applied the zero- 
energy resolution where needed. Exclusion limits calcu- 
lated with omission of the energy-dependent terms give 
more conservative limits than would be obtained with in- 
clusion of those terms, and omission allows the Si and Ge 
detectors to be treated in the same manner. 

As will be discussed in more detail in Section IIV CI 
we made use of two forms of zero-energy resolution when 
calculating WIMP exclusion limits: the YNR-corrected 
and the Q-corrected recoil-energy resolutions. The for- 
mer is needed to include the effect of subthreshold pulses 
occasionally triggering the experiment due to phonon- 
noise fluctuations, and therefore does not include any 
noise contributions from the ionization channel. The lat- 
ter represents the resolution intrinsic to the combined 
efficiency of the analysis cuts, and includes ionization as 
well as phonon noise. The other-detector triggered data 
were used to measure the YNR-corrected and Q-corrected 
recoil-energy resolutions at keV. The results are listed 
in Table UTT1 for the viable low-threshold detectors. 

Though the electron-recoil calibration is valuable for 
measuring detector performance, the nuclear-recoil en- 
ergy scale is most important to WIMP searches. We 
characterized this response by comparing the spectrum 
of nuclear recoils observed in 252 Cf calibrations with a 
GEANT3 [27j simulation of the source and detector geom- 
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FIG. 7. (color online). Comparison of Monte Carlo simulated 
recoil-energy spectra (blue/dark error bars) to experimen- 
tally measured recoil-energy spectra (orange/light error bars) 
from 252 Cf calibrations for the Ge (top panel) and Si (bottom 
panel) detector ensembles. Exponential fits are overlaid for 
both data (blue/dark solid) and Monte Carlo (orange/light 
solid). Despite the clear disagreement in absolute event rate, 
comparison of the fitted decay constants provides a calibra- 
tion of the nuclear-recoil energy scale for each target mate- 
rial. There is a 16±3% discrepancy in the Si energy scale. A 
smaller discrepancy is observed in the Ge energy scale, and is 
within two standard deviations (statistical) of expectation. 



etry. In order to improve statistical power, we summed 
the observed events to produce mean Ge and Si recoil 
spectra. Figure [7J shows the Monte Carlo and data dis- 
tributions for both detector types with exponential fits, 
where the Si (Ge) spectra are well described by one (two) 
decaying exponential(s). Comparison of the fitted decay 
constants tests our nuclear-recoil energy scale without 
being sensitive to the obvious disagreement in absolute 
rate. Although difficult to perceive in Fig. [Jj there is a 
substantial 16±3% discrepancy in the Si detectors' en- 
ergy scale. Without distinctive energy peaks with which 
to confirm the energy scale, we resolved the discrepancy 
by adjusting the Si recoil-energy scale upwards by 16%, 
which has the effect of making the Si upper limit (dis- 
cussed in Section HVCj) more conservative. The Si hard- 
ware and software thresholds and analysis cut efficiencies 
discussed above reflect this correction. A smaller discrep- 
ancy exists for the Ge detectors, where the low-energy 
portion of the recoil distribution is overestimated relative 
to the Monte Carlo, and the high-energy portion is un- 
derestimated. We did not make an analogous correction 
because the discrepancy is smaller (within two standard 
deviations of expectation) , and the calibration done with 
electron-recoil peaks is more reliable. The choice to not 
correct the Ge recoil-energy scale leads to a slightly more 
conservative exclusion limit for low-mass WIMPs. 

The effect "ion channeling" might have on a direct 
detection experiment's energy scale has received con- 
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TABLE IV. The number of raw live days and correspond- 
ing number of WIMP candidate events for each detector and 
WIMP search following application of the analysis cuts and 
energy thresholds is listed. The smaller exposure and corre- 
sponding number of candidates for the Z2 6V data is due to 
the event burst cut. 



3V Exposure 6V Exposure 

Detector Live Days Candidates Live Days Candidates 



Z2 


66.12 


159 


20.16 


67 


Z3 


66.12 


129 


51.66 


349 


Zl 


66.12 


130 


51.66 


125 


Z5 


66.12 


174 


51.66 


202 


Z6 


66.12 


401 


51.66 


314 




Recoil energy (keV) 



siderable attention within the dark matter community 
since Sekiya et al.'s [28| 2003 publication. Early mod- 
els [29|, [3(| for a Nal scintillation detector suggested that 
a potentially significant fraction of recoiling nuclei might 
be "channeled" between rows and planes within a detec- 
tor's crystalline lattice, allowing them to interact more 
weakly with atomic cores and deposit more energy in 
the form of ionization. The full recoil energy from a 
WIMP interaction might therefore be detected via elec- 
tron recoils, beneficially lowering a Nal detector's effec- 
tive energy threshold for nuclear recoils. In our detec- 
tors, the larger ionization signal associated with chan- 
neling would cause an elevated ionization yield for nu- 
clear recoils, tending to push channeled nuclear recoils 
closer to the electron-recoil population. Ion channeling 
would therefore have a negative effect upon our nuclear- 
recoil acceptance efficiency. If the effect is particul arly 
pronounced for low recoil energies as suggested in [30| . 
the low-mass WIMP sensitivity of this analysis could be 
severely impacted. 

An upper limit on the fraction of channeled events near 
threshold can be calculated with the 252 Cf data plotted 
in Fig. 2J For recoil energies between 2 keV and 6 keV, 
the fraction of events with ionization yield above the 
nuclear- recoil band is less than 4%. Although this rate is 
higher than the rate of electron recoils observed without 
the source present, neutron interactions with materials 
near the detectors can lead to secondary gamma rays. 
In particular, neutron captures on hydrogen in the in- 
ner polyethylene shielding yield a continuum of electron- 
recoil energies due to the 2.2 MeV photons released in the 
process, biasing our channeling upper limit high. There 
were also several ~lMeV photons emitted directly from 
the 252 Cf source with each fission. A limit on the chan- 
neling fraction that is less conservative, by incorporat- 
ing a Monte Carlo estimate of the rate of nonchanneled 
electron recoils during 252 Cf calibrations, is beyond the 
scope of this paper. Nevertheless, we do not see signif- 
icant evidence of a channeling effect in our data that is 



FIG. 8. (color online). The combined Ge (top panel) and 
Si (bottom panel) WIMP candidate event rates as a function 
of recoil energy. The uncorrected event rates (blue/dark thin 
solid) are compared to the efficiency corrected event rates. 
The latter are successively corrected by the exposure- weighted 
detector-averaged efficiencies of the analysis cuts (black/dark 
thick solid), and then by the detector-averaged hardware and 
software thresholds (orange/light thin solid). 



large enough to appreciably affect the efficiencies esti- 
mated in this analysis. We therefore ignore the effect of 
ion channeling, a decision that is supported by the recent 
and more sophisticated models developed by Bozorgnia 
et al. (3~lj ]. which indicate that ion channeling for cryo- 
genic Ge and Si targets is effectively nonexistent for low 
recoil energies. 



IV. RESULTS 

A. Candidate events 

Following application of the analysis cuts and phonon 
software thresholds, a substantial residual rate of events 
is observed in the low-threshold signal region. We re- 
strict our attention to events with recoil energies between 
0.5 keV (IkeV) and 100 keV for the Ge (Si) detectors, 
yielding a total of 1080 Ge and 970 Si candidate events. 
The number of raw live days and corresponding number 
of candidate events for each detector and bias voltage are 
listed in Table ITVl The combined Ge and Si recoil-energy 
spectra are shown in Fig. [8l where the event rates have 
been successively corrected by the average efficiencies for 
the analysis cuts, and then by the average hardware and 
software thresholds. Since the recoil-energy spectra and 
the former efficiencies are functions of Q-corrected recoil 
energy, before dividing out the latter efficiencies they are 
converted from iNR-corrected to Q-corrected recoil en- 
ergy by smearing with the ionization noise. 
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B. Backgrounds 

Although the recoil spectra resemble in shape the dis- 
tributions expected for WIMP interactions, the events in 
the signal region are likely due to several types of unre- 
lated background processes consisting of electron recoils, 
zero-ionization events, 14 C contamination particular to 
Z6, and nuclear recoils from cosmogenic neutrons. We 
will not subtract these events, but will accept them as 
candidates for the purpose of calculating upper limits on 
a WIMP signal; this is the most conservative treatment 
of these data. 

A few background populations are particularly evident 
in plots of ionization yield versus recoil energy. The sig- 
nal regions and candidate events for representative Ge 
and Si detectors are displayed in Fig. [9] The most easily 
identified background is specific to the Ge detectors. The 
distinct 1.3 keV line between 1 keV and 3 keV in recoil en- 
ergy accounts for a substantial number of the candidate 
events. On average, the internal electron capture x-rays 
or Auger-electrons from the decays of 68 Ge and 71 Ge have 
unit ionization yield. Due to the relatively low signal-to- 
noise ratio at these energies in both phonons and ion- 
ization, however, electronic noise induces a large tail of 
electron recoils to low ionization yield. The feature is 
tilted with respect to the recoil-energy axis because of 
anticorrelation between the numerator and denominator 
of the ionization yield expression. Leaked electron recoils 
from the 1.3 keV line account for ^20% of the Ge can- 
didates in the 3V WIMP search, and for approximately 
one-third (Z3) to one-half (Z2 and Z5) of the Ge candi- 
dates in the 6V WIMP search. The fraction is greater for 
the 6V data because extensive 252 Cf neutron calibrations 
were performed prior to this data period, enhancing the 
levels of 71 Ge via thermal neutron capture. 

Although more prominent for the Si detectors, both 
the Ge and Si detectors are afflicted by a distribution 
of events with nearly zero ionization yield. Close exam- 
ination reveals that these "zero-charge" events possess 
ionization signals indistinguishable from electronic noise. 
They are otherwise normal events with recoil energies 
as large as 100 keV in some instances. The distribution 
of zero-charge events for the Si detector in Fig. [9] spans 
the full analysis energy range, and clearly crosses into 
the signal region for recoil energies <7keV, constitut- 
ing the majority of this detector's candidate events. A 
possible explanation is that these events were the result 
of recoils that occurred near the side edges of the ZIP 
detectors, where electric field lines did not span the de- 
tectors' 1 cm thicknesses. Ionization occurring along field 
lines that terminate on a detector's edge rather than on 
an electrode is not properly drifted across the crystal, 
resulting in no signals in either the inner or outer ion- 
ization electrodes. For recoil energies >10keV, the zero- 
charge events' xy-positions can be reliably reconstructed, 
and tend to cluster near detector edges. Past analy- 
ses of CDMS data have avoided zero-charge events by 
requiring WIMP candidates to have a minimum ioniza- 




Recoil energy (keV) 

FIG. 9. (color online). Ionization yield versus recoil energy for 
unvetoed single scatters passing the data-quality and fiducial- 
volume cuts (all dots) for representative Ge (top panel, Z5 
6V) and Si (bottom panel, Z4 3V) WIMP searches. Each de- 
tector's signal region is outlined in this plane by its nuclear- 
recoil band (blue/dark solid), phonon energy software thresh- 
old (yellow/light solid), and the extent of the horizontal axis. 
Although these regions are partially cut off from above and 
below, all 202 Z5 6V and 130 Z4 3V WIMP candidate events 
(red/dark dots) can be seen. 



tion energy, analogous to the phonon software thresholds 
described above. Unfortunately, an ionization threshold 
would severely limit our WIMP detection efficiency for 
recoil energies <5 keV. 

The significantly higher event rate for Z6 is believed to 
be due to 14 C surface contamination. Prior to Run 21, 
during testing at one of the CDMS test facilities, Z6 was 
operated in close proximity to a detector that had been 
previously exposed to a 14 C calibration source with faulty 
encapsulation, accidentally contaminating one of its sur- 
faces with a low level of the isotope. For this reason, 
Z6 was placed at the bottom of the detector tower with 
its contaminated surface facing away from the adjacent 
detector. Beta decays of 14 C produce electrons with an 
average energy of ~50keV and a maximum energy of 
^156 keV. Beta radiation in this energy range will in- 
teract entirely within a ZIP detector's ~10/xm surface 
dead layer, where the charge collection efficiency is con- 
siderably reduced. These events have reduced ionization 
yield, and populate the gap between the bands of elec- 
tron and nuclear recoils when plotted in the fashion of 
Fig. [9l with a substantial number leaking into the signal 
region. Although surface events can be rejected with high 
efficiency for recoil energies >10keV through a combina- 
tion of phonon and ionization pulse timing parameters, 
the near-threshold WIMP detection efficiency cannot be 
preserved. 

Leakage of electron recoils into the nuclear-recoil band 
is a component of each detector's candidate events, al- 
though the source is usually Compton scatters of pho- 
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tons. The discrimination based on ionization yield 
breaks down as the recoil energy decreases, until at the 
"crossover energy" the electron- and nuclear-recoil bands 
significantly overlap. The crossover energy varies from 
~3keV (Z5 3V) to -7keV (Z6 6V). More Neganov- 
Luke phonons are produced when the detectors are run 
at higher bias voltage, causing degradation in both the 
recoil-energy resolution and the yield-based discrimina- 
tion. Consequently, ZIP detectors perform better as low- 
threshold detectors with the lower 3V bias voltage. Scal- 
ing calibration data from a 60 Co source results in the 
estimate that only a few Compton electron-recoil events 
per detector leak into the nuclear-recoil band for recoil 
energies above the crossover energy. We have not de- 
vised a reliable method of estimating the contribution of 
Compton electron-recoil leakage for recoil energies below 
the crossover energy. We estimate that 10% to 20% of 
the WIMP candidates are actually electron recoils from 
Compton scatters. 

The highest-energy signal events are largely due to the 
neutron background associated with the SUF's modest 
overburden. Muons (and hadronic showers produced by 
them) occasionally broke apart nuclei in the rock sur- 
rounding the experiment, expelling high-energy neutrons 
with sufficient energy to punch through our shielding 
and create lower-energy neutron secondaries within the 
shielding materials capable of producing signals above 
threshold. The expected rate of these unvetoed neu- 
tron interactions as a function of recoil energy was sim- 
ulated for the Ge detectors (dashed curve in Fig. 2 in 
18]). With no efficiencies applied, the rate peaks at 
~0.1 events kg _1 keV _1 d _1 at the 0.5 keV cutoff, and de- 
cays quasiexponentially to ~0. 003 events kg _1 keV _1 d _1 
at 100 keV. Application of the Ge detectors' average de- 
tection efficiencies, followed by a scaling to the ~71 kg- 
days of Ge exposure, yields an expectation of ~66 neu- 
trons among the 1080 Ge candidate events. 

An accurate accounting of the contribution of each 
source of background to the total number of candidate 
events is difficult. Most candidate events have low recoil 
energies for which our experimental variables are unable 
to differentiate among the various sources, particularly 
for the zero-charge and Compton electron-recoil leakage 
backgrounds. An approximate tally of the percentage of 
candidate events due to each background source is listed 
in Table [V] for the combined Ge and Si detector ensem- 
bles, where an "other" category is included to indicate 
the percentage of events that are not attributed to the 
background sources as described above. The other events 
could be due to unidentified sources of background or 
WIMPs, but probably reflect systematic uncertainty in 
the estimates for the known background sources. 



C. Exclusion limits 

The large uncertainties associated with our back- 
ground sources preclude the subtraction of backgrounds. 



TABLE V. Percentages of the 1080 Ge and 970 Si WIMP 
candidate events caused by each of the known background 
sources: l)leaked electron recoils from the 1.3 keV line; 2)zero- 
charge events occurring near detector edges; 3)leaked electron 
recoils due to beta decays of 14 C embedded in the surface of 
the Si detector Z6; 4) leaked electron recoils due to Compton 
scatters of photons; and 5)nuclear recoils due to neutrons. 
The final category labeled "other" indicates the percentage 
of events not attributed to the known background sources as 
estimated. 



Ge (%) 



Si (%) 



Electron Capture 1.3 keV Line 32 

Zero-Charge Events 30-40 30-40 
14 C Contamination Betas 40 

Compton Photon Electron Recoils 10-20 10-20 
Cosmogenically-Induced Neutrons 6 2 

Other 2-22 0-18 



We can claim no evidence of a WIMP signal, and pro- 
ceed to calculate conservative limits under the assump- 
tion that the candidate events may constitute a WIMP 
signal. We use the observed event rates to set upper lim- 
its on the WIMP-nucleon spin-independent elastic scat- 
tering cross section as a function of WIMP mass. Two 
limits are calculated; the main result of this paper is a 
combined Ge and Si exclusion limit based upon all 2050 
signal events, while a secondary result focuses on the Si 
data alone. 

To calculate exclusion limits the observed event rates 
must be compared to a hypothetical WIMP model. For 
convenient comparison with the results of other direct 
detection experiments, we work within the framework of 
the "standard" halo model described in [3] , but normal- 
ized to a local WIMP density of 0.3GeV/cm 3 . We also 
conformed to the dark matter community's standard as- 
sumptions for the WIMP characteristic and mean Earth 
velocities of 220km/s and 232km/s, respectively [llij . 
while assuming the more recently estimated value of 
544km/s for the galactic escape velocity jlEj . The effect 
of the galactic escape velocity on these results is discussed 
in more detail in Section IIV Dl 

Limits were calculated using a version of Yellin's "op- 
timum interval" method [HJ that has been extended to 
accommodate high statistics [33j]. For each choice of 
WIMP mass, the limit is essentially determined by a 
single energy interval for which the number of observed 
events is particularly low relative to the expected num- 
ber of WIMP events. An appropriate statistical penalty 
is applied for choosing the interval that sets the best 
limit, yielding a 90% confidence level upper limit. This 
method is especially effective at discriminating against 
backgrounds which are distributed differently from the 
expected signal. 
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To calculate a single exclusion limit, the data from the 
individual Ge and Si detectors have to be appropriately 
combined. Traditionally, CDMS has combined the detec- 
tor ensemble into a single averaged detector, where the 
individual-detector masses and efficiencies are averaged 
according to their exposures. This "averaged" method 
for combining detectors makes use of the entire exposure, 
and comingles the candidate event energies for different 
detectors before forming the energy intervals required by 
the optimum interval method. Figure |S] represents the 
averaged version of the candidate event data for this anal- 
ysis. The averaging method is appropriate when the de- 
tectors involved have approximately equal sensitivity to 
WIMP interactions, as was the case for previous CDMS 
WIMP-search results in which the analyses were either 
background-free or nearly so. 

When the averaging technique is applied to detectors 
with variable event rates, the detectors with especially 
high event rates effectively pollute the lower-rate detec- 
tors by filling in the most sensitive intervals with a dispro- 
portionate number of events. For this reason, we decided 
to adopt a novel "serialization" technique for combining 
the detector data. Energy intervals are separately pre- 
pared for each detector in order to preserve the most 
sensitive intervals. The intervals are then concatenated 
in an arbitrary order which, to avoid possible bias, was 
selected before the effect of the order was known. We 
chose to place the 3V data before the 6V data, and then 
to order them according to their position within the de- 
tector tower (from top to bottom). If the limit-setting 
intervals do not span multiple detectors, the order will 
not affect the result. This technique allows the opti- 
mum interval method to calculate the limit from the best 
individual-detector energy intervals. The resulting limit 
reflects only a fraction of the exposure, rather than the 
total exposure for the entire detector ensemble. This is 
a trade-off we decided to accept before calculating the 
limits. Each detector is clearly background-limited, par- 
ticularly near threshold where our low-mass WIMP sensi- 
tivity resides. Trading exposure for cleaner energy inter- 
vals should yield stronger limits for low masses. In hind- 
sight, this turned out to be true for WIMP masses less 
than 8GeV/c 2 . The serialization technique also allows 
for different detector types within the detector ensemble, 
providing a natural method for combining the Ge and Si 
data. 

To include the effect of non-zero energy resolution 
properly, expected WIMP rates were separately calcu- 
lated for each detector and WIMP search (3V and 6V 
data) in a series of steps. The limit was calculated for 
75 WIMP masses between 1 GeV/c 2 and 100GeV/c 2 . 
At each mass, the halo model predicts the differen- 
tial WIMP-nucleon scattering rate in terms of an ideal, 
perfect-resolution recoil energy (see Fig. Q] for example). 
Each detector's ideal spectrum was then convolved with 
its YNR-corrected recoil-energy resolution listed in Ta- 
ble Mil (first two columns). Recall that the hardware 
trigger and software phonon thresholds depend solely 
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FIG. 10. (color online). Comparison of 90% confi- 
dence level upper limits from the combined Ge and Si 
(black/dark solid, our main result) and Si only (gray/light 
solid) data, with potential signal regions based on data from 
the DAMA/LIBRA [H and CoGeNT [H experiments. The 
two (larger) oval-shaped filled regions (pink/light shaded) 
represent the DAMA/LIBRA annular modulation signal as 
interpreted by Savage et al. (99.7% C.L.), and include 
the effect of ion channeling as modeled by Bozorgnia et 
al. [3l]]. The (smaller) oval-shaped filled region is the 
99% (gray /medium shaded) confidence level signal region 
found by Hooper et al.'s [36| simultaneous best fit to the 
DAMA/LIBRA and CoGeNT data. The elongated filled re- 
gions are SUSY theory predictions by Bottino et al. [l| 
for J7wimp < ficDMmin (dark-yellow/medium shaded) and 
£2wimp > ficDMmin (blue/dark shaded). Our limits assume a 
galactic escape velocity of 544km/s [lSj. while the potential 
signal regions are based on a value of 600km/s. 



on the phonon signal. The expected WIMP spectrum 
should therefore include noise from only the phonon 
channel when the hardware and software threshold ef- 
ficiencies are applied. After application of these phonon- 
only efficiencies, the spectrum was further smeared to 
include the electronic noise of the ionization channel via 
a second convolution with the quadrature difference be- 
tween the Q-corrected (third and fourth columns of Ta- 
ble [nlj and YNR-corrected recoil-energy resolutions. The 
threshold-reduced expected WIMP spectrum, in terms 
of Q-corrected recoil energy as measured by a ZIP detec- 
tor, was then multiplied by the remaining analysis cut 
efficiencies, which either depend weakly on this energy 
estimator or are constant. Finally each detector's dou- 
bly smeared and efficiency-reduced expected WIMP rate 
was scaled by the detector's mass and exposure (listed 
in Tables U and IIV[) . The resulting distribution of ex- 
pected events per Q-corrected keV for each detector was 
used to construct a cumulative probability that describes 
how likely it is that a WIMP interaction deposited an 
energy within each energy interval defined by that de- 
tector's candidate event energies. The 10 (4) individ- 
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FIG. 11. (color online). Comparison of 90% confidence 
level upper limits from the combined Ge and Si (black/dark 
solid, our main result) and Si only (gray /light solid) data, 
with those from XENON100 with constant (orange/light 
dotted) or decreasing (orange/light dashed) scintillation ef- 
ficiency extrapolations at low energy [33, CoGeNT [3^] 
(+), and CRESST [H (blue/dark dash-dotted). Our lim- 
its (and XENONlOO's) assume a galactic escape velocity of 
544 km/s OH , while the CRESST limit uses 650 km/s and the 
CoGeNT limit uses 600 km/s. See also Fig.[T2]in which limits 
for other escape velocities are compared. 



ual probability distributions were serialized as described 
above, and a combined Ge and Si (Si only) optimum in- 
terval upper limit was calculated at each WIMP mass. 
Note that although particular care was taken to include 
the effect of non-zero energy resolution for this result, it 
would have made only a trivial contribution to the sen- 
sitivities of previous CDMS results due to their higher 
thresholds. 

The resulting exclusion limits are shown in Fig. I1Q[ 
and partially exclude parameter space associated with 
potential signals from the DAM A/LIBRA [U, HI and 
CoGeNT [35( experiments. The former includes the 
(small) effect of ion channeling as modeled by Bozorgnia 
et al. [3l| . The combined Ge and Si limit cuts through 
the middle of the 99% confidence level signal region as- 
sociated with Hooper et al. 's 36] simultaneous fit to the 
DAM A/LIBRA and CoGeNT data, and excludes new pa- 
rameter space for WIMP masses between 3GeV/c 2 and 
4GeV/c 2 . Our limits are compared to those from other 
experiments with strong low-mass WIMP sensitivity in 
Fig. 1111 The additional information required to repro- 
duce our results can be found in [3{|, and includes the 
efficiencies and candidate event energies for each detector 
and voltage bias. 

A limit based on the Ge data alone was also calculated, 
and found to be almost equal to the combined Ge and 
Si limit (the difference is imperceptible in Fig.s [TU] and 
[T2")) . The serialization technique sets the limit according 



to the most sensitive detectors. The Ge detectors were 
more sensitive across the range of WIMP masses consid- 
ered, yielding a combined Ge and Si limit that is almost 
the same as the Ge only limit. We decided beforehand 
that the combined detector limit would be our main re- 
sult, but that we would publish the Ge only and Si only 
limits as well. The limit for low masses (<6GeV/c 2 ) 
was determined by Z2 3V and 6V intervals, while a com- 
bination of the Z5 3V, Z5 6V, and Z2 3V data set the 
limit for WIMP masses >10GcV/c 2 . The sensitivity in 
the intervening mass range was dominated by the Z3 3V 
data. Each Si detector contributes to the Si only sen- 
sitivity, with the Z6 data providing the best limit for 
masses <4GeV/c 2 , and the Z4 data for larger WIMP 
masses. None of the intervals chosen by this serialized 
implementation of the optimum interval method spanned 
more than one detector, confirming that the order of the 
serialization is unimportant. 



D. Systematic studies 

A recent study by Smith et al. [15j analyzed the high- 
est velocity halo stars in an early release of data from 
the RAVE survey [4(j, and derived a loose but convinc- 
ing constraint on the galactic escape velocity. At the 
90% confidence level, their preferred range of escape ve- 
locities is 498 km/s to 608 km/s, with a median value 
of 544 km/s. Though based on an otherwise identically 
normalized standard halo model, the results for the other 
experiments shown in Figs . 1 1 Ol and 1111 assume a variety of 
galactic escape velocities which are generally larger. The 
CRESST [38[ limit is based on a much larger 650 km/s 
value, while the DAM A/LIBRA and CoGeNT results (in- 
cluding the Hooper et al. region) use 600 km/s. Although 
we have chosen to adopt the more recent measurement of 
544 km/s for this paper, additional limits were calculated 
to see the effect of the larger, 650 km/s escape velocity. 
There is very little visible difference in the main result 
for WIMP masses larger than 4 GeV/c 2 , while the Si only 
limit assuming the lower escape velocity is as much as 
20% weaker for WIMP masses between 6GeV/c 2 and 
10GeV/c 2 . The effect of changing the escape velocity 
becomes important for masses less than 4GeV/c 2 . Lim- 
its based on the 650 km/s and 544 km/s galactic escape 
velocities for WIMP masses down to lGeV/c 2 are com- 
pared in Figure 1121 Regions corresponding to Smith et 
al.'s 90% confidence interval are plotted as well. The ef- 
fect is nearly maximal for a 2 GeV/c 2 WIMP mass, where 
the combined Ge and Si limits span over an order of mag- 
nitude, and the Si only limits span over two orders of 
magnitude. 

The detection efficiencies associated with each detec- 
tor's hardware and software phonon thresholds, and anal- 
ysis cuts, were estimated directly from data. In each case, 
the estimates are based on a finite number of events, 
resulting in statistical uncertainties. These uncertain- 
ties have a small (but noticeable) effect on our experi- 
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FIG. 12. (color online). The 90% C.L. upper limits from 
the combined Ge and Si and Si only data for WIMP masses 
down to 1 GeV/c 2 and an assumed galactic escape velocity of 
544km/s [l5| (black/dark solid and gray/light solid, respec- 
tively), compared to limits from the same data for an escape 
velocity of 650km/s 16] (black/dark dashed and gray/light 
dashed, respectively). The ranges of combined Ge and Si 
(blue/dark shaded) and Si only (yellow/light shaded) limits 
corresponding to the 90% C.L. range of escape velocities found 
by Smith et al. 15!], and the limits from the CRESST [H 
(blue/dark dash-dotted, 650km/s) and CoGeNT [H (+, 
600km/s) experiments are shown as well. 



mental sensitivity for WIMP masses less than 6GeV/c 2 . 
The combined Ge and Si exclusion limit for a 5GeV/c 2 
WIMP, for example, is -20% weaker when the 90% C.L. 
lower limit efficiencies (1.28ct less than the mean) arc 
substituted for the best-fit mean efficiencies. The ef- 
fect is maximal for a 1 GeV/c 2 WIMP, where the same 
WIMP exclusion limit is a factor of 5 weaker. Again, 
we decided to adopt a conservative approach and em- 
ploy the 90% C.L. (statistical) lower limit efficiencies in 
all our WIMP exclusion limit calculations. The best-fit 
mean efficiencies are compared to these reduced efficien- 
cies in Fig. [5] for representative Ge and Si detectors. All 
figures and tables in this paper preceding Figs. [10] [12] re- 
flect the best-fit mean efficiencies where applicable, while 
Figs. [TUHT21 are based on reduced detection efficiencies 
such as those shown in Fig. [3] 

The combined Ge and Si exclusion limit was also 
tested for sensitivity to uncertainties in the Ge detec- 
tors' electron-recoil energy scales. Each Ge ZIP's Ier- 
corrected and Q- correc ted recoil energy scales for elec- 
tron recoils were checked for the entire 3V and 6V WIMP 
searches using the lines described previously. The three 
peaks' mean energies were estimated with Gaussian fits 
relative to energy scales originally derived from short 
137 Cs calibrations. Deviations between the observed 
and expected peak energies were generally within sta- 
tistical fluctuations. A few detectors, however, exhib- 
ited sizable phonon signal overestimates for the 1.3 keV 



line. The worst case was for the Z2 6V data, where the 
YER-corrected recoil energy was over measured by —6%, 
and the Q-corrected recoil energy was over measured by 
—16%. The latter is worse because it includes an ion- 
ization signal under-measurement in conjunction with a 
phonon signal over- measurement. An ad hoc recalibra- 
tion of the Ge recoil-energy estimators was performed 
and the WIMP exclusion limits were recalculated. Be- 
cause the energy scale discrepancies mostly involve recoil- 
energy overestimates, the combined Ge and Si limit based 
on the recalibrated energy scales is slightly stronger. The 
limit is mostly insensitive to the recalibration, however, 
with the difference in limits only becoming visible by eye 
for the lowest WIMP masses. The effect is maximal for 
a WIMP mass of 1 GeV/c 2 , where the recalibrated limit 
is a factor of —2 stronger. The final limits shown in 
Figs. [T0H121 are based on the original Ge detectors' en- 
ergy scales, a slightly conservative choice that avoids the 
systematic uncertainty associated with the ad hoc recal- 
ibration. 



It is easy to imagine that our low-mass WIMP sensitiv- 
ity might be heavily reliant on the lowest-energy trigger 
efficiency, and on our inclusion of the effect of non-zero 
energy resolution near threshold. If true, then relatively 
modest unknown systematic uncertainties in the hard- 
ware trigger efficiencies could lead to large uncertain- 
ties in the WIMP exclusion limits. The combined Ge 
and Si limit was tested for sensitivity to potential trig- 
ger efficiency systematic uncertainty in two ways. First, 
each detector's hardware trigger efficiency was set equal 
to zero wherever it drops below 10%, eliminating the 
lowest-energy detection efficiency. The limit based on 
the reduced trigger efficiencies is reassuringly not much 
weaker than the limit using the full efficiencies. The 
difference is maximal at lGeV/c 2 , where the reduced- 
efficiency limit is —3% weaker. Second, the sensitivity to 
possible systematic uncertainty associated with the in- 
clusion of near-threshold energy resolution was tested by 
smearing the hardware trigger efficiencies by 50% of the 
^NR-corrected recoil-energy resolutions (~50eV for Ge, 
and — lOOeV for Si) before applying them in the limit 
calculation described above. This addresses a concern 
that there might be an unknown extra resolution factor 
included in the hardware trigger efficiency estimates due 
to the difference between the phonon signal that triggers 
the experiment and the offline reconstructed phonon en- 
ergy. It would actually be more appropriate to decon- 
volve the trigger efficiencies, but the reverse operation is 
easier to implement, and although it will give a stronger 
rather than weaker limit, the size of the effect should be 
similar. The resulting limit is indeed stronger, but only 
slightly so. The size of the effect reaches a maximum for 
a lGeV/c 2 WIMP mass, where the difference is —5%. 
Our low-mass WIMP sensitivity appears to be relatively 
insensitive to systematic uncertainties in the hardware 
trigger efficiencies. 
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V. CONCLUSIONS 

Within the dark matter community, possible signals 
from the DAM A/LIBRA and CoGeNT experiments have 
aroused renewed interest in relatively low-mass WIMPs. 
Despite substantial backgrounds, CDMS ZIPs are ex- 
cellent low-threshold detectors. The analysis presented 
here is background-limited and conservative, yet sets 
competitive limits for low WIMP masses on the spin- 
independent WIMP-nucleon elastic scattering cross sec- 
tion. Under standard halo assumptions, we partially ex- 
clude the parameter space favored by interpretations of 
the DAM A/LIBRA annual modulation and CoGeNT ex- 
cess as WIMP signals. Ignoring the effect due to uncer- 
tainty related to the galactic escape velocity, the limits 
are robust for WIMP masses greater than ~2GeV/c 2 . 

The methods developed for this low-threshold analy- 
sis are not limited to the Run 21 shallow-site data. An 
upcoming analysis based on CDMS II data taken at our 
deep site in the Soudan Mine employs a similar but lower- 
background analysis technique, and achieves similar 
recoil-energy thresholds for a significantly greater expo- 
sure. We have also experimented with operating CDMS 
detectors at very high bias voltages (~100 V), sacrificing 
our two-channel discrimination in favor of ultra-low en- 
ergy thresholds (~50eV) via Neganov-Luke phonon am- 
plification [4lJ. The CDMS II data are a rich source 



of information for probing low-mass WIMP parameter 
space. 
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